Magneto-transport of large CVD-grown graphene 
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We present magnetoresistance measurements on large scale monolayer graphene grown by chemical 
vapor deposition (CVD) on copper. The graphene layer was transferred onto SiCh/Si via PMMA 
and thermal release tape for transport measurements. The resulting centimeter-sized graphene 
samples were measured at temperatures down to 30mK in a magnetic field. We observe a very 
sharp peak in resistance at zero field, which is well fitted by weak localization theory as well as 
strong localization. The samples exhibit conductance fluctuations symmetric in field, which are 
due to large scale inhomogeneities consistent with the grain boundaries of copper during the CVD 
growth. 

PACS numbers: 81.05.ue, 73.22.Pr, 72.80.Vp, 81.15.Gh, 73.20.Fz 



Graphene has attracted a considerable amount of at- 
tention due to the ease in isolating a single sheet of 
graphite via mechanical exfoliation [l], Q. Despite the 
fact that it is one atom thick, exfoliated graphene has 
shown extraordinary transport properties and can be 
used as a novel material for many potential applications. 
Its unique band structure has led to many interesting 
phenomena such as tunable charge carriers between elec- 
trons and holes, anomalous integer quantum Hall effect 
P, IH and ultrahigh mobilities at room temperature [1] . 
The most noteworthy property of the band structure is 
the existence of two degenerate Dirac cones [H, which 
leads to two degenerate valleys (K and K'). 

In a seminal work by McCann and co-workers @ based 
on earlier work on honeycomb lattices Q, the authors 
have obtained a general expression for the weak local- 
ization (WL) correction in graphene, which determines 
the dependence of the magnetoresistivity as a function 
of the magnetic field (B) for various scattering param- 
eters in graphene. The increase in resistance at zero 
field is due to an increased backscattering probability 
when a path and its time reversed path is phase coherent. 
This time reversal symmetry is destroyed by a small B, 
which leads to a decrease in resistance with field. This 
generic WL effect is suppressed by inelastic scattering 
due to a loss of phase coherence. In graphene, however, 
backscattering can only occur when scattering between 
the two valleys is possible (inter- valley scattering). This 
depends on the type of scattering potential. Typically, 
slow varying potentials will lead to no inter- valley scatter- 
ing, whereas short range impurities will. The existence of 
inter-valley scattering and intra- valley scattering (within 
a valley) leads to an interesting competition between WL 
and weak anti-localization (WAL), since in the absence 
of inter-valley scattering, only WAL exists, and is re- 
sponsible for the increase in resistance with B. Intention- 
ally disordered exfoliated graphene promotes inter- valle 
scattering, which can then lead to strong localization 

While early experiments on exfoliated graphene have 
shown a strong suppression of WL even at very low tem- 



peratures Q, more recent experiments have shown sig- 
nificant WL effects by averaging over many carrier con- 
13- Hj]. This was necessary because meso- 



centrations 

scopic conductance fluctuations dominate transport in 
these small graphene flakes since the sample sizes are 
comparable to the phase coherence length ~ 1/im 
at low temperatures. Transport in large scale graphene 
would be self-averaging when the sample size is much 
larger than Lq. This is what we propose to demonstrate 
in this letter. 

While acquiring large graphene flakes by exfoliation 
might not be feasible, other techniques have been de- 
veloped for large-scale production, such as epitaxial 
graphene grown on a SiC(OOOl) surface, where WAL has 
been observed [l3| . Recently, nickel films [lij and copper 
foils |l5[ have been successfully used as substrates to pro- 
duce large scale single layer graphene by chemical vapor 
deposition (CVD). 

For this work we grew graphene monolayers by CVD 
of hydrocarbons on 25 /xm-thick commercial Cu foils. 
The Cu foil is first acid-treated for 10 mins using acetic 
acid and then washed thoroughly with de-ionized water. 
Graphene growth is realized in conditions similar to Li 
et al. 15 1 . but using a vertical quartz tube. Graphene is 
grown at 1025°C in 0.5 Torr, with a 4 seem H2 flow and 
a 40 seem CH4 flow for 30 minutes. The methane flow 
is stopped while the hydrogen flow is kept on during the 
cooling process. 

The graphene film was transferred onto n ++ -doped Si 
wafer with an oxide thickness of 285nm. This is realized 
by etching the Cu foil in an oxidizing solution of 0.1M 
ammonium persulfate ((NH4) 2 S 2 8 ). PMMA or Nitto 
Denko thermal release tape (3195MS) are used as a sup- 
porting layer. For some devices, we did not remove the 
PMMA layer unlike Li et al. [l(| but used it to provide an 
additional insulation between the gate and the graphene 
as well as to keep the integrity of the film. As for the ther- 
mal release tape, its adhesive force is removed by heat- 
ing the tape to 120°C in order to detach the graphene 
film and deposit it onto Si02/Si. Small contacts of SPI 
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FIG. 1: Left: transfer process of the graphene monolayer. Top 
right: SEM picture of graphene on copper showing typical 
grain boundaries. Bottom right: picture of the lx 0.5cm 2 
graphene device on a SiCb/Si substrate used in this work 
with mm-sized silver paste contacts. 



conductive silver paint, arranged according to the van 
der Pauw technique, are placed on the periphery of the 
graphene film and on the Si backgate. Silver wires are 
used to connect the contacts to the chip carrier to per- 
form transport measurements. An image of the sample 
is shown in fig. [T] Dozen cm 2 -sized samples were made 
and all were conducting, with four terminal resistivities 
ranging between 1 and 4kf2 at room temperature. All 
samples show a weak increase in resistance upon cooling 
to 4.2K. Contact resistances vary between 1 and 10kf2. 
Raman spectroscopy and optical reflection were used to 
confirm the monolayer nature of graphene on the devices 



17]. 



The samples were placed inside a dilution refrigerator 
and measured as a function of B. A most striking feature 
is a sharp peak of the resistance at B = as shown in fig. 
[21 Similar peaks were obtained in all samples measured. 

Because we have no lithographically defined contacts, 
which avoids possible residues from lithography, there 
is a large symmetric contribution in the Hall trace due 
to contact misalignment. However, this large symmetric 
signal stemming from the relatively large magnetoresis- 
tance (~2kf2) can be effectively eliminated by antisym- 
metrization. The result is a remarkably linear antisym- 
metric Hall trace (see fig. which shows no indication 
of quantum Hall plateaus. The slope can be used in order 
to effectively extract the total Hall carrier density, i.e., 
n = {edR^/dB)- 1 ~ 7.3 • 10 12 cm- 2 . By applying a 
gate voltage we observe a change in the Hall slope (den- 
sity) of the sample. In the remainder we will use a zero 
gate voltage configuration for a more detailed analysis. 



FIG. 2: Top graph: Magnetoresistance (R xx ) as a func- 
tion of B at 30mK. Also shown are the extracted symmet- 
ric (R x v x m = \R XX {B) + R xx (-B)]/2) and antisymmetric 
(R xx ym = [R XX (B) - R xx (-B)]/2) traces. Bottom graph: 
Hall resistance (R xy ) as a function of B and the corresponding 
symmetric and antisymmetric components. The Hall slope is 
i?^ m =-(85.29 ± 0.005)fi/T. 



This (anti-)symmctrization method described above can 
also be applied to the magnetoresistance trace in fig. [2l 
where we observe a very small linear antisymmetric con- 
tribution. This small antisymmetric Hall contribution in 
R xx comes from a non-ideal contact region. In most cases 
this contribution can be neglected. 

We now turn to understanding the large and very nar- 
row peak at B = 0, which strongly suggests WL. Follow- 
ing the theoretical work by McCann et al. we can 
write the magnetoresistivity as a function of B as: 



p(B)=p(0)— P —{F{B^)- 



-F( J B +2S i )-2F(S + J B i +-B*)), 

(1) 

where F(z) = ^(0.5 + z/B) - \n{z/B), 4> is the 
digamma function and p = h/e 2 is the quantum unit 
of resistance. B a with a = (j>, i, or * is the char- 
acteristic B associated with the scattering time T a — 
ey/imp(0)/2TTVFB ( j, where vf — 10 6 m/s is the graphene 
Fermi velocity, is the dephasing time, tj the inter- 
valley scattering time and r* the intra-valley scattering 
time. We can now use equation (1) in order to fit our 
experimental data. The result is shown in fig. [3l 

Overall, the quality of the fit is quite remarkable con- 
sidering that the data is obtained without perfor ming any 
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averaging unlike the case for exfoliated graphene 10h12 | . 
While £>0 is strongly determined by the curvature very 
close to B = 0, the fitting parameters Bi and B* are 
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FIG. 3: In the bottom graph we show the symmetric re- 
sistance as a function of B for different temperatures (30- 
600mK). The WL fit is shown with a dashed line and offset 
by +10n for clarity and without the B ~ region. The top 
left graph shows the same data plus temperatures from 1.1 
to 23K but zoomed in around zero field. The top right graph 
shows the temperature dependence of the results to the best 
fits from equation (1) for parameters B,p, Bi, B+ and i?(0) 
with lines as guides for the eye. The lowest two temperatures 
were measured one week apart from the other temperatures 
and show very small changes in the pattern of the resistance 
fluctuations. 
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TABLE I: Values of various parameters for the lowest tem- 
perature data (30mK) for two values of the resistivity. The 
mobilities p are in cm 2 /V-s and the scattering lengths L CT in 
nm. The geometric factor of the sample gives a resistivity of 
p=(3±0.5)kf2 and the fits are good over a wider range (up to 
6T) than for a resistivity of lkfi (up to 4T). 95% confidence 
intervals for fits of the scattering times are similar for the two 
resistivities and are shown in figure [3] 



less well defined and strongly depend on the range of 
fields over which the fit is performed. Therefore, we first 
perform a fit over a wide range in B (-6 to 6T), which 
strongly determines B* and then fit for and Bi over 
a much narrower range (-.08 to .08T) using B+ obtained 
from the first fit. The temperature dependence of the 
parameters is shown in fig. [3j An important additional 
parameter in the theory is the value of the resistivity at 
zero field p(0). In particular, n is strongly dependent on 
the value of p(0). We illustrate this in table 1, where 
we show the values of various parameters for p(0) = lkf2 
and p(0) = 3kQ. 



The mobility is given by p = l/enp(0), the trans- 
port time by r tr = pa/2p{Q)^JwnvF @, the mean free 
path by L tr = vpT tr and the other scattering lengths by 
L a = VFy / T a T tr /2 (<t = cj), £, or *). While the physical di- 
mension of the contact configuration on the sample leads 
to a resistivity of ~ 3kO, it could well be that the actual 
relevant resistivity is significantly smaller if the sample is 
inhomogeneous. For instance, p(0)=lkf2 corresponds to 
the lower bound from the different samples we measured. 
In the presence of large inhomogencities, the current will 
meander along the current path with the least resistance 
and hence the actual current path might be significantly 
longer than the direct one, which would reduce the ef- 
fective resistivity compared to the geometrically defined 
one. Hence the values of Bi and Ti have to be taken with 
some caution as they are very dependent on the value of 
resistivity. For instance, a change of factor 3 in resistivity 
can affect Tj by more than a factor 10 because the best 
fits for Bi give a difference of a factor 20. Interestingly, 
when only considering scattering lengths, differences are 
not as marked (factor 2.5) as shown in table 1. 

In general, Bi determines the field at which the nega- 
tive magnetoresistance due to WL stops and turns over to 
WAL. If Bi = there is no WL and only WAL is visible, 
whereas for Bi —¥ oo only WL survives. Hence the ratio 
Bi/B^ sensitively determines the height of the WL peak. 
It is interesting to compare our scattering lengths with 
those obtained for exfoliated graphene. For instance, ex- 
periments on exfoliated graphene have obtained values 
for Li ranging between 75 and 250nm [10Hl2j . whereas 
for the maximum these experiments obtained values 
ranging between 0.6 and lfim and for values between 
20 and lOOnm. Overall, the values are quite similar and of 
the same order, with > Lj > L+ at the lowest temper- 
atures. With increasing temperature, only decreases 
due to increased inelastic scattering, while the other scat- 
tering lengths remain largely temperature independent. 
There is also a monotonous decrease in resistance with 
increasing temperature, which cannot be explained by 
WL alone, since WL would predict a much smaller effect 
1- 

It is therefore instructive to analyze the resistance peak 
at B = in terms of strong localization. Indeed, the 
localization length, L c , in graphene is determined by a 
"universal" B-dependence (independent on disorder and 
size) of AL c ((f))/ L c ((p) as shown in the inset of figure |U 
AL c (0) = L c (4>) - L c (0) and = B ■ voia{Ll,L\} is 
the amount of flux within the area determined by L? c or 
Li. For <C L c and assuming a localized behavior 
p r^, e L 4>/ L " ~ 1 + L^/Lc, the dip in the localization 
length at B = is given by AL C /L C (B) ~ -0.4 (see 
fig- SI), which leads to a peak in the resistivity given by 
Ap/p(0) ~ O.4L /L C (O). The 30mK data can then be 
nicely fitted with this expression and yields L^/L c (0) ~ 
0.11 (see figure [3}, hence L c (0) ~ ~ 4pm, when 
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FIG. 4: (a) Inset: AL C /L C as a function of the magnetic 
flux. The different data points correspond to different values 
of disorder and sample widths (W=40-160 in units of carbon 
atoms), which all collapse on the same curve. Main figure: 
Ap/p as a function of B for the 30mK data (straight line) 
and for p ~ exp(a/L c ) (points) with a/L c (B — 0) ~ 0.11. 
The error bars correspond to the standard error of the mean 
obtained from the numerical calculation of L c for graphene 
ribbons of widths W. (b) Resistance fluctuations as a function 
of the magnetic field shown as the difference from the WL fit. 



using = 450nm from table 1. This could explain the 
observed increase in resistance beyond the WL prediction 
due to the proximity of strong localization. can also 
be estimated directly from the B ■ L^-dependence of AL C , 
noting that the width of the dip is determined by BL 2 ^ ~ 
h/e). This yields ~ 550nm, which is consistent with 
L<j, obtained using McCann's theory. 

We now turn to the fluctuations seen in figures El and 
SJ In mesoscopic samples (L ~ L$) such as in exfoliated 
graphene, mesoscopic conductance fluctuations dominate 
transport properties [ll|. Here we have a macroscopic 
sample where the sample is much larger than the coher- 
ence length, which averages out mesoscopic conductance 
fluctuations. So the remaining fluctuations are of a dif- 
ferent nature, i.e., classical. Classical fluctuations can 
arise from particular disorder distributions, like ensemble 
averaged conductance fluctuations [HI or as more likely 
here, from large scale inhomogeneities. We can estimate 
a lower bound on the size of the inhomogeneous regions, 
by assuming that the inhomogeneities are uncorrelated. 



For fluctuations of the order of 1%, as seen in figures [3] 
and[H this would lead to L corr > 0.01L, where L ~ 1cm 
hence L corr > 100/mi. Interestingly, the typical grain size 
of the Cu surface used for the CVD growth is of the order 
of 500/iin (see figure [TJ, which we measured using elec- 
tron backscatter diffraction (EBSD). Since the growth of 
graphene is affected by the crystal orientation, this could 
lead to large scale inhomogeneities with graphene grown 
on different grains having different densities and mobil- 
ities, which would explain the observed fluctuations of 
the resistance in B. These fluctuations are perfectly sym- 
metric in B due to the small Hall resistance and are re- 
producible for a given contact configuration even after 
several weeks. However, they change drastically when 
changing the density by only 0.5%, when applying a small 
voltage on the gate formed from the doped Si substrate as 
seen in fig. |4j This also implies that it would be difficult 
to observe quantum oscillations at low fields. Interest- 
ingly, these large inhomogeneities can be circumvented 
by patterning contacts on small parts of the sample as 
was done in ref . [n| , where the quantum Hall effect was 
observed. 

Summarizing, we have studied the magnetotransport 
of large scale graphene monolayers obtained by CVD 
growth and subsequent transfer, which show many in- 
teresting features. Most importantly, they exhibit a very 
sharp peak at B = in resistance, which can be at- 
tributed to WL as well as strong localization. We fur- 
ther observed large resistance fluctuations, which might 
be due to the macroscopic grain boundaries from the orig- 
inal CVD growth and demonstrated an effective way to 
extract very precise Hall densities even in the presence 
of these large resistance fluctuations and large contact 
misalignments. 
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